On the large-scale instability in interacting dark energy and dark matter fluids 
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Recently, Valiviita et al. (2008) have reported a large-scale early-time instability in coupled 
dark energy and dark matter models. We take the same form of energy-momentum exchange and 
specialise to the case when the interaction rate is proportional to Hubble's parameter and the dark 
energy density only. Provided the coupling is made small enough for a given equation of state 
parameter, we show that the instability can be avoided. Expressions are derived for non-adiabatic 
modes on super-horizon scales in both the radiation and matter dominated regimes. We also examine 
the growth of dark matter perturbations in the sub- horizon limit. There we find that the coupling 
has almost no effect upon the growth of structure before dark energy begins to dominate. Once the 
universe begins to accelerate, the relative dark matter density fluctuations not only cease to grow 
as in uncoupled models, but actually decay as the universe continues to expand. 



I. INTRODUCTION 

There is good evidence to believe the present day en- 
ergy density of the universe is mostly in the form of dark 
energy [3, HI, the properties of which remain relatively un- 
known. Furthermore, observations of Type la supernovae 
[3, H| leave little doubt that the expansion of the universe 
is accelerating. Viable models of cosmology now require 
a large dark energy component, capable of producing the 
negative pressure required for accelerated expansion. 

By far the simplest model of dark energy is Einstein's 
cosmological constant, A. The cosmological constant(A) 
and cold dark matter (CDM) model, with values of to- 
day's density parameter for the dark energy U, x w 0.7 and 
dark matter f2 m w 0.3 is the current prevailing paradigm. 
But while consistent with observational constraints, the 
standard model is in many ways unsatisfactory. One such 
example is the 'coincidence problem': why are the energy 
densities in the dark energy and dark matter comparable 
today, when the redshift dependence of each is so differ- 
ent? 

Motivated to explain the coincidence problem while 
deviating as little as possible from the successful ACDM 
model, a coupling between dark energy and dark matter 
has often been considered. An energy exchange modifies 
the background evolution of the dark sector, and explain- 
ing the coincidence problem can be reduced to tuning a 
coupling parameter to an appropriate value. 

The coupling enters via the continuity equations. With 
energy exchange rate Q between the dark energy (sub- 
script x) and the cold dark matter (subscript c), the dark 
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energy obeys the continuity equation in conformal time 

p x + 3H(l+w x )p x = -Q, (1) 

while the dark matter obeys 

p c + m Pc = Q. (2) 

Here we have introduced the equation of state parameter 
wa that gives the ratio of the pressure Pa to the energy 
density pa of a fluid, 

w A = — • (3) 

PA 

We have also used H = aH, where a(t) is the expan- 
sion scale-factor and H the Hubble parameter. Acceler- 
ation of the expansion rate requires the energy density 
of the universe to dominated by a fluid with an effective 
equation of state parameter w c g < —1/3. We do not 
allow the phantom case of w < — 1 in this work. Sim- 
ple solutions for the background exist for couplings of 
the form Q = aTLp x + (3Tlp c . These were initially in- 
vestigated by Chimento [|[ and then expanded upon by 
Barrow and Clifton who provided general solutions 
for any cosmology with two components exchanging en- 
ergy in such a fashion, provided the components were 
modelled as cosmological fluids with constant w. Quar- 
tin et al. [3] examined the observational constraints upon 
such a class of models, significantly limiting the available 
parameter space. Again, the equation of state parameter 
was treated as fixed. Non-zero values of (3 were found to 
reduce the required fine-tuning of the initial energy den- 
sity, as well as increase the observationally allowed values 
of a 0]. 

A coupling would influence more than just the back- 
ground dynamics of the universe. In particular, the 
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growth of perturbations in the coupled fluids would be af- 
fected. Recent work by Valiviita et al. Q has shown that 
couplings of the simple form described above, with con- 
stant w, exhibit extremely rapid growth of dark energy 
fluctuations on super-horizon scales in the early universe. 
In fact, the perturbations in the dark energy become un- 
stable for any model with non-zero /3, no matter how 
small this parameter is made. While this would appear 
to rule out all couplings of the above form and with con- 
stant w, the explicit examples in [H included no cases 
where the interaction rate was proportional to the den- 
sity of dark energy and not of the dark matter, i.e. with 
(3 — and a ^ 0. Here we look at just such a scenario. 



II. BACKGROUND EVOLUTION 

Friedmann's equation relates the evolution of the scale- 
factor a(t) to the background energy density p. We make 
use of conformal time, r, which is related to cosmic time 
via dt = adr. Overdots indicate derivatives with respect 
to conformal time. 

Friedmann's equation reads 



8ttG 



-pa . 



(4) 



With the choice of Q = aTtp x , the continuity equations 
can be solved to yield [a, la] 



Px — Px,0 a 



-(3(l+™)+a) 



^ = aPxfl a _( 3 (l +t0 )+a) 



3w + a 



Pc,0 



aPx,o 
3w + a 



-3 



(5) 



(6) 



We follow the standard notation where a subscript zero 
indicates today's value. We normalise the scale-factor so 
that do = 1. The ratio of dark energy to dark matter 
density r can then be written 



Pr_ 

Px 



Pc,0 
Px,0 



3w 



3w+a 



3w + a 



(7) 



With 1 3it; | < a, the dark energy and dark matter ap- 
proach a constant ratio as the universe expands. The 
coincidence problem can be said to be solved if this ratio 
is of order unity, but this requires a value of a already 
observationally excluded Q- Nevertheless, as argued in 
, non-zero values of a can still be said to alleviate the 
problem. We restrict ourselves to positive values of a. 



III. PERTURBED FRW COSMOLOGY 

We assume a flat FRW cosmology and work in Newto- 
nian gauge, 

- ds 2 = dt 2 (l + 2*) -a 2 (l- 2^)S ij dx i dx j , (8) 



with metric signature (+,—,—,—). We work in Fourier 
space, using comoving Fourier wave- vectors k % = hi, so 
that d l d i -> k 2 /a 2 . 

The four-velocity of fluid A is given by 



u w = {0--*)> a ~ lv U 



(9) 



The peculiar velocity three-vector v l = w, are small. We 
define the velocity perturbation 6 = diV 1 . 



A. Energy-momentum tensors 

The energy-momentum tensor for fluid A is given by: 
T »tA) = (/A) + p(A)\ i;(A) MU (A) _ 5 ^ P (A)^ (10) 

The total energy-momentum tensor is simply the sum of 
the components, 



(11) 



We define the density perturbation in fluid 8a using 
Pa = (1 + 5a)Pa- An overbar denotes the background 
quantity, though we will usually leave this implicit. 

Energy and momentum conservation for fluid A implies 



(A) 



and conservation for the entire system requires 



(12) 



E^)=°- ( 13 ) 



The four-vector Q^a) governs the energy exchange be- 
tween components, and it is to this we now turn our 
attention. 



B. Covariant energy exchange 

The energy exchange in the background does not de- 
termine a fully covariant form of energy exchange 0, 0] . 
Instead, an energy exchange four-vector must be spec- 
ified. We adopt the approach of Q and consider two 
scenarios; aligning the four-vector with the dark energy 
four-velocity, 



Q(A) ~ QaU£ , 



(14) 



or with the four- vector of the dark matter four- velocity, 



Q{A) - 



(15) 



These choices produce slightly different outcomes, and 
the differences are noted as we proceed. 

To produce the desired changes to the continuity equa- 
tions, we see that aQ c = —aQ x — aHp x in both cases. 
We also make the common assumption that aH gives 
an interaction rate that has no spatial dependence. We 
therefore perturb only p Xl not 7i, in the coupling. 
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C. Sound speed of dark energy 

The speed of sound of a fluid or scalar field A is denoted 
by c s a- For a barotropic fluid with a constant value of 
u>a, then c 2 A = wa- This leads to an imaginary speed of 
sound for the dark energy (c 2 x = w x < 0) . An imaginary 
sound speed leads to instabilities in the dark energy; the 
problem is commonly remedied by imposing a real sound 
speed by hand. A common choice (and the one we make 
here) is the scalar field value of c sx = 1. 

This choice leads to an intrinsic non-adiabatic pressure 
perturbation in the dark energy. This contains a term, 
highlighted recently in Q , that arises due to the coupling 
between dark energy and dark matter. We include this 
term, and refer the interested reader to Q . 



D. Perturbation equations of motion 

Conservation of the energy-momentum tensor, com- 
bined with results of the previous sections and our choice 
of energy exchange four-vector, implies the following. For 
the dark energy density perturbation: 



s' x + m{i-w x )5 x + (i + w x )e x + m 2 {i-w 2 x )^ 



3(l + w x )& = -aH 



tf + 3W(l-0-§ 



(16) 



For the dark energy velocity perturbation, the right-hand 
side differs slightly depending on our choice of energy 
exchange four-vector. 



0' - 2H9 X 



1 + to. 
(l + b)aH 
l + w x 



8 X - kH 



where 



6 = 



if = QaU? 

1 ifQf A) = 



X 1 



(17) 



(18) 



For the dark matter, the density perturbation obeys 
S' c + e o - 3$' = oM P — [6 X - S c ] , (19) 

Pc 

while the velocity perturbation is governed by 

e' c + e c n-kH 

= (i-b)aH^[9 x -9 c }. . (20) 

Pc 

The perturbed Einstein equations are well known, and 
we do not reproduce them here. They can be found in 
[Toj . whose notation for the scalar metric perturbations 
we share. 



IV. INITIAL CONDITIONS IN THE EARLY 
RADIATION ERA 

In Q it was shown that models with (3 ^ suffered 
from an early time large-scale instability no matter how 
small the value of (3. This was driven by a term propor- 
tional to (3 on the right-hand side of equation (jTTJ) . A 
term proportional to a also exists, which can be large if 
w is close to —1 or a is made very large. In this section 
we examine how large this term needs to be to cause the 
non-adiabatic mode to be a growing one. 

We consider super-horizon scales (k/Ti. <§; 1) and as- 
sume adiabatic initial conditions. The gravitational po- 
tentials are dominated by fluctuations in the dominant 
fluid (radiation or matter). The well known result is that 
$ oc ^ — constant. The constant of proportionality in 
the radiation era is determined by the anisotropic stress 
generated by the neutrinos. In the absence of neutrinos 
or in the matter dominated era, the potentials are equal. 
These assumptions will be invalid only if perturbations 
in the dark energy are large enough to influence the grav- 
itational potentials. As the dark energy has a very low 
background density in the radiation era, this can only 
happen if 5 X grows extremely large. 

Neglecting time derivatives of the gravitational poten- 
tial, and keeping only leading order terms in k/H, the 
dark energy equations p^|) - (jTTJ) can be combined into 
a second order equation: 



n i - 3w 



(l + b)a 



1 



= (AH 2 + BH') (21) 

The constants A and B have values unimportant for our 
analysis. 

In the radiation era, 7i = t -1 . The adiabatic mode is 
therefore an obvious solution: S x oc = constant. To 
find the remaining solutions, we define a new variable 
S x = $x + C^P, with the constant C chosen such that the 
right-hand side of (|2"Tj) is equal to zero. In the radiation 
dominated era, we can then write: 



t 2 5" + 3 - 3m; 



(l + 6)a 
1 + w 



+ 3(1 - w) 2 - b 



l + w 



5 X = (22) 



When 6=1, equation (|2"2"|) becomes formally the same 
equation found by He et al.[ll|, despite the differing as- 
sumptions made about the physics involved. In their in- 
vestigation of perturbations given a background coupling 
of the form Q = aTLp x , they choose to set the net mo- 
mentum exchange to zero (Q^a) = 0), m contrast to our 
adoption of the form of momentum exchange used in [f| . 
The differences between the 6=1 choice of momentum 
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exchange and zero net momentum exchange arise in the 
equations for the dark matter perturbations, which are 
not used in the above analysis, nor in the analysis by He 
et al. This leads to the same behaviour of dark energy 
perturbations. This is not true when 6 = 0, and can re- 
sult in different behaviour (oscillatory or non-oscillatory) 
for the same choice of parameters (see the remainder of 
this section). Note also that the simplifying assumptions, 
and their justifications, made in 11] differ to those made 
here: we have neglected terms that will be small due to 
choice of intial conditions, and simplified the result by 
extracting the adiabatic mode. In [ll| , terms are instead 
neglected that are found to be small from a numerical 
analysis. 

are power laws, S x oc r n± . 



Solutions of equation J22 
The index is given by: 



n± = 



± 



2(1 + to) 2(1 + io)' 



(23) 



where we follow the notation of [ll| and have defined the 
quantities 



and 



3w 2 + to + (1 + b)a - 2, 



A = 9to 4 + 30w 3 + (13 - 6(6 - l)a)w 2 
+ 2w [(1 + b)a - 14] + 4(26 - l)a 
+ (l + 6) 2 a 2 -20. 



(24) 



(25) 



In the limit of w very close to -1 (and assuming a is 
reasonably small), we can expand as a series in (1 + to), 



2(1 + to) 



-5/2- 



(1 + b)a 3(1 + to) 
2(1 + to) + 2 

+ 0(l + w) 2 , (26) 




FIG. 1: The evolution of the dark energy outside the horizon 
in a matter dominated universe, for modes of three different 
scales. We take a = 0.08 and w = —0.98, and 6 = 1. The 
agreement with the analytical approximation is excellent until 
the mode begins to leave the horizon (kr ~ 1). Vertical lines 
indicate when this occurs for each mode. 



and 



A « 4(1 + 36)a 2 + 12(56 - 3)a(l + to) 

+ (24(1 - 6)a-62)(l + w) 2 + 0(1 + w) 3 . (29) 

Once again, the second term in T can result in a rapidly 
growing dark energy fluctuation. 

We have solved equations ([T6"]) - (|20p numerically in 
the matter dominated regime (Figure [1]) , where we need 
not worry about the radiation fluid and its perturba- 
tions. The analytical agreement is excellent until the 
mode leaves the horizon (kr ~ 1). Numerically we sec 
that when this happens the mode begins to oscillate with 
a growing amplitude. 



A w (l + 36)a 2 + 2(76-5)(l + w)a 

+ (6(l-6)a-23)(l + w) 2 + C(l+w) 3 . (27) 

When a = 0, the non-adiabatic mode is decaying. But 
when the coupling is switched on, the second term in 
r can become very large, resulting in n + 3> 1. For a 
range of a and to, which is much larger in the 6 = case, 
oscillatory behaviour can also result (due to A becoming 
negative). The instability means these coupled models 
suffer from all the problems outlined in @] for /3 / 
models, unless the value of a is made small enough. The 
closer to is to —1, the smaller a must be made to avoid the 
instability. This is in contrast to /3 ^ models, which are 
unstable no matter how small the parameter /3 is made. 

In the matter dominated era, we can carry out the 
same procedure, this time with Ti. = 2r~ 1 . We find, 

on F x ~ -9/2+^±-^+3(l+ W )+0(l+ w ) 2 , (28) 



V. SUB-HORIZON EVOLUTION IN THE 
MATTER AND RADIATION DOMINATED ERAS 

In the sub- horizon limit, H 2 /k 2 <C 1, equation (fT6|) 
yields, 

6' x + 3W(1 - w x )5 x + (1 + w x )6 x - 3(1 + w x )& = 0. (30) 

Note the two terms on the right-hand side of equation 
(fTtj)) scale as H 2 /k 2 . As these are the only two terms con- 
taining the coupling parameter a, the simplified equation 
above does not contain the coupling parameter. 

One of the perturbed Einstein equations simplifies to 
Poisson's equation in comoving coordinates, 

- fc 2 * = 4^Ga 2 (p x 6 x + Pc S c ) . (31) 

Without the coupling, the dark energy perturbations are 
significantly suppressed on small-scales in comparison to 
dark matter perturbations, primarily due to its large 
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speed of sound [12j |. The coupling does nothing to alter 
this fact unless the right-hand side of equation (fT?)) makes 
a significant contribution. If the early time instability has 
been avoided this cannot be the case, as a/(l+w) will be 
small. Thus it is reasonable to expect the dark energy to 
remain suppressed on sub-horizon scales. We therefore 
neglect dark energy perturbations for the remainder of 
this section. 



By combining equations (|19p and (|20[) , we eliminate 9 C 
and find a second-order equation for the growth of the 
matter density perturbation. From the above argument, 
we have neglected dark energy perturbations. Keeping 
only the dominant gravitational terms, 



2a^ ) 5' 



n ( i 

Pc 

a— (H' - U 2 (a + 3w - 1)) 5 C = -fc 2 *. (32) 

Pc 



We note that in the limit of a — > 0, this reduces to the 
standard growth equation, with the well known growing 
mode S c oc t 2 in both matter and radiation eras. The ad- 
ditional terms are proportional to ar (recall r is the ratio 
of dark energy to dark matter). In the matter dominated 
regime, then ar <C 1, and these terms will be negligible. 
Even when r *~ 1, the terms will be suppressed by the 
size of a, which will be small itself. The dominant effect 
causing a deviation from standard linear growth of struc- 
ture in the matter dominated regime will therefore be, 
as in an uncoupled cosmology, the influence of the dark 
energy upon the expansion rate. The growth of struc- 
ture in a coupled model can therefore be treated in the 
matter dominated regime simply as an uncoupled model 
with an effective dark energy equation of state parame- 
ter w e s = w + a/3. This will cease to be true only when 
the background energy density of matter is no longer well 
approximated by its usual p c oc a~ 3 dependence, and the 
late time scaling behaviour becomes apparent. 



with the new time variable r\ = r — T m . Note that as 
rj increases (r decreases and approaches Too)) the scale- 
factor increases. The constant of integration, t^, is the 
radius of the de Sitter event-horizon in the uncoupled 
case with a cosmological constant. The growth equation 
can then be written as, 



2 trl, 



2 - 12w - 4a 
a + 3w + 1 



6', 



^ (3a + 9w — l)(3w + a) + aw 

(a + 3w + l) 2 c 



0. (34) 



This admits power law solutions, 8 C oc r) m , where 




8a 



w(l + 3w + a) 5 



(35) 



In the range of a and w relevant to the problem, then 
m > 0. Recalling that r\ decreases with increasing scale- 
factor, we see that the universe becomes steadily more 
homogeneous as it expands. We interpret this to be a 
combination of two effects. The first is the accelerating 
expansion, which slows and (without the coupling) even- 
tually stops structure formation. This occurs, for exam- 
ple, in ACDM cosmology when the cosmological constant 
becomes dominant. The second effect is that dark en- 
ergy is constantly being transformed into dark matter, 
via the coupling. As the rate is proportional to the den- 
sity of the dark energy, and the dark energy density is 
essentially uniform, new dark matter is also created uni- 
formally. This rising 'background' of dark matter reduces 
the relative value of the fluctuations, reducing 5 C . 

We have also investigated both numerically and ana- 
lytically the extreme late time behaviour, where krj -C 1 
and the modes can be thought of as leaving the horizon. 
We find the tend toward homogeneity continues, but with 
a much milder rate of decay. 



VII. CONCLUSIONS 



VI. SUB-HORIZON EVOLUTION IN THE 
DARK ENERGY DOMINATED ERA 



The coupling between dark energy and dark matter 
eventually leads to a constant ratio between the two dark 
components. With a small value of a, the dark energy 
still dominates. We consider the evolution of structure 
once this equilibrium has been reached. 



Friedmann's equation solves to yield 



H = 2(a + 3w + l)~ 1 r)- 



(33) 



We have shown that constant w models with the same 
form of energy-momentum exchange considered by Q 
suffer from an instability with a ^= 0, even if /3 = 0. 
However the instabilities in these models are not as se- 
vere as those facing models with (3^0. There is at least 
some non-trivial region of parameter space where the in- 
stability can be avoided, although the value of a is now 
constrained both from background observables Q and 
from stability requirements to be extremely small. De- 
spite this, any non-zero value of a will lead to a late-time 
scaling regime, alleviating (even if not solving) the coinci- 
dence problem. It is unfortunate that with a constrained 
to such small values, we find any observable trace upon 
the growth of CDM structure will be negligible. Detect- 
ing a coupling of this form from measurements of large- 
scale structure is extremely doubtful, even with the pre- 
cision promised by future experiments. 
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We have said nothing up to this point of models of dark 
energy with a variable equation of state parameter, such 
as scalar-field (quintessence) models. The same caveats 
in 0] apply here. Much of the above analysis will not 
apply in variable w models, although some parameteri- 
sations such as the often used w = wq + (1 — a)w a lead to 
fixed w over large periods of time. Our analysis will apply 
during those epochs of constant w. We refer interested 
readers to recent work on quintessence with couplings of 
this or similar form (such as recent work [l3l Il4j). 

The future decay of dark matter fluctuations is an in- 
teresting result. It implies observers today find them- 
selves close to the time of maximum inhomogeneity. The 
more the coincidence problem is alleviated, the closer to 
the late time scaling regime today becomes, and thus the 
closer to the peak of inhomogeneity. Without the cou- 
pling, observers find themselves at the time of the end 
of structure growth. The root cause in both cases is the 



acceleration of the universe only beginning today. Our 
position as apparently privileged observers in this fashion 
remains difficult to explain in any satisfactory way. 
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Note added in proof: after submission we became aware 
of the work by [1 51 ] in which an independant analysis of 
general instabilities in coupled models was carried out, 
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for negative values of a and find a non-trivial region con- 
sistent with observations. 
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